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ABSTRACT. Dynamic degrees of freedom and internal variables are treated in a uniform 
way. The unification is achieved by means of the introduction of a dual internal variable. 
This duality provides the corresponding evolution equations depending on whether the 
Onsager-Casimir reciprocity relations are satisfied or not. 
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1. Introduction 

Introduction of internal variables suggests taking into account the influence of an inter- 
nal structure on the dynamic behaviour of a material fl~). As long as evolution equations of 
internal variables are determined by general macroscopic principles, we can expect that the 
validity of the evolution equations is independent of particular microscopic models. It is 
important, therefore, to understand the universal aspects behind these evolution equations. 

There are two basic methods to generate the evolution equations for internal variables. 
Both methods are based on fundamental principles. 

The first method generates the evolution equations exploiting the entropy inequality. 
This approach uses exclusively thermodynamic principles, and the corresponding variables 
are called internal variables of state Q. This frame has the advantage of operating with 
familiar thermodynamic concepts (thermodynamic force, entropy), however, no inertial 
effects are considered. The thermodynamic theory of internal variables has a rich history 
(see the historical notes in |3|). A first more or less complete thermodynamic theory was 
suggested by Coleman and Gurtin ||4], and the clear presentation of the general ideas of the 
theory was given by Muschik [5 1. Internal variables were applied for several phenomena 
in different areas of physics, biology, and material sciences. A complete description of 
the thermodynamic theory with plenty of applications based on this concept of internal 
variables of state can be found in @. 

The second method generates the kinetic relations through the Hamiltonian variational 
principle and suggests that inertial effects are unavoidable. This approach has a mechanical 
flavor, and the corresponding variables are called internal degrees of freedom. Dissipation 
is added by dissipation potentials. This theoretical frame has the advantage of operating 
with familiar mechanical concepts (force, energy). The method was suggested by Maugin 
0, and it also has a large number of applications |8]|9). The clear distinction between 
these two methods with a number of application areas is given by Maugin and Muschik 
EHU1 and Maugin 0. 

It is important to remark that the terminology of the field is not unique, and even more, 
it is contradictory. For example, Verhas [6] calls his internal variables as dynamic degrees 
of freedom following the suggestion of Kluitenberg (see e.g. |[TTl[T2l ). Sometimes internal 
variables or internal degrees of freedom appear without thermodynamic foundation under 
different names. An important example can be found in nerve pulse dynamics, where the 
classical "phenomenological variables" of Hodgkin and Huxley [13] and the "recovery 
variables" of Nagumo, Arimoto and Yoshizawa lfl4l [l5 1 are internal variables from a ther- 
modynamic point of view, as it was shown in 1 16 1. Another important example is damage 
mechanics (e.g. [ 17|), where the macroscopic damage variable is an internal variable from 
a thermodynamic point of view, and the thermodynamic framework can reveal several new 
properties of the theory fl8l Q1D . 

In this paper we follow the terminology of Maugin and Muschik 1 3 1 with some impor- 
tant extensions. We call internal variables of state those physical field quantities - beyond 
the classical ones - whose evolution is determined by thermodynamical principles. We call 
internal degrees of freedom those physical quantities - beyond the classical ones - whose 
dynamics is determined by mechanical principles. 

One of the questions concerning this doubled theoretical frame is related to common 
application of variational principles and thermodynamics. Basic physical equations of 
thermodynamical origin do not have variational formulations, at least without any further 
ado 1201 . That is well reflected by the appearance of dissipation potentials as separate 
theoretical entities in variational models dealing with dissipation. 
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On the other hand, with pure thermodynamical methods - in the internal variables ap- 
proach - inertial effects are not considered. Therefore, the coupling to simplest mechanical 
processes seemingly requires to introduce some improvements, those are usually new prin- 
ciples of mechanical origin. 

Therefore, one can have the impression that the doubling of the theoretical structure is 
a necessity, because the usage of both mechanical and thermodynamical principles can- 
not be avoided. This is the conceptual standpoint of the GENERIC approach (GEneral 
Non-equilibrium Equation of Reversible Irreversible Coupling) [21 1. However, a doubled 
theoretical structure is not very convenient, because it doubles the number of physical as- 
sumptions restricting predictive capabilities of the theory. 

In what follows, we propose a uniform approach based exclusively on thermodynamic 
principles. Our suggestion requires dual internal variables and a generalization of the usual 
postulates of non-equilibrium thermodynamics: we do not require the satisfaction of the 
Onsagerian reciprocity relations. With dual internal variables we are able to include in- 
ertial effects and to reproduce the evolution of dynamic degrees of freedom. It could be 
impossible with a single internal variable. This is the price we pay for the generalization. 
In other words, instead of the doubling of the theoretical structure we suggest the doubling 
of the number of internal variables. 

In the next Section we give a short overview of the thermodynamic frame, then we 
summarize the theory of internal variables and internal degrees of freedom in Sections 3 
and 4, respectively. The unified frame is developed in Section 5. Finally, we shortly discuss 
the advantages and disadvantages of the suggested approach. 



2. Thermodynamics of continua 

First, we recall the basic equations of the finite-strain thermoviscoelastic continua with 
internal variables. In this case, the state space of a material point X is spanned by the 
variables (e, F, a), where e is the specific internal energy, F is the deformation gradient 
(direct-motion gradient) defined as F = J^- | t , and a. = (a±, a n ) are internal 

variables of state. Here x = x(X, t) denotes the mapping from the reference configuration 
to an actual one. 

The balance of linear momentum in an actual configuration reads (in the absence of 
body forces) 11221 

(1) pv-Vt = 0, 

where v = % | x is the velocity field, the dot above the velocity denotes the substantial 
time derivative, and t is the Cauchy stress. The density p is not independent of the deforma- 
tion gradient, because the density in the reference configuration is po = p dct F = const. 

The balance of internal energy can be calculated as the difference between the conserved 
total energy and the kinetic energy and is given as (e.g. (6)) 

(2) pe + V-q = t : FF 1 . 

Here q is the heat flux and FF 1 = Vv. Then the entropy balance can be represented in 
the form (6] El 

(3) ps + V • j = q • Vi + - (t - t™») : FF 1 - 1a • d. 



4 



P. VAN*, A. BEREZOVSKI AND J. ENGELBRECHT 



Here we introduced the following traditional thermoelastic quantities (given by Gibbs's 
equation): 

(ds ds ds 

-p—,-p—,...,-p— 
oai oot2 oa n 

The entropy s as a function of the state space variables (e, F,a) is a concave function 
according to thermodynamic stability requirements. 

It is important to emphasize that the entropy flux j in Eq. ([3]) is assumed to have its 
classical form 

(6) 3d = j;- 

This is a severe assumption, and in the case of mixtures and weakly nonlocal extension of 
classical theories the entropy flux deviates from this classical form as j = y + J. This 
additional term was proposed first by Miiller ll23ll . The particular form of the additional 
term J depends on the thermodynamic interactions and their level of nonlocality (e.g., in 
the case of weakly nonlocal heat conduction a family of its different forms is presented in 
El). 

In what follows we remind the derivation of the evolution equations of internal variables 
from general principles. As the coupling to mechanical or thermal interactions could be 
involved and would blur the structure of equations, in the following we restrict ourselves 
to the case, where the internal variables are completely decoupled from the mechanical 
and thermal interactions. Namely, this is the situation where the mechanical interaction is 
reversible t = t rev , the heat flux is zero q = 0, p = const., v = const., and there is 
a complete decoupling at the static level as well, that is s(e, F, a) = st{c) + sjvf(F) + 
sa{oi). Let us remark that one can give very different conditions of decoupling, especially 
if other thermodynamic potentials are introduced, as, e.g., the Helmholtz free energy. 

3. Internal variables of state 

Let the basic thermodynamic state be spanned by n scalar internal variables a = 
(ai, a.2, a n ). We want to determine their evolution according to thermodynamic prin- 
ciples. 

3.1. Local state. Let us assume that the evolution equation of a. is given in a general form 

(7) a = g a , 

where the dot denotes the material time derivative that can be regarded as a partial one in 
our investigations. 

The local state hypothesis of Kestin Il25ll26ll is introduced by the assumption that en- 
tropy and the right hand side of the evolution equation (0 are functions of a and do not 
depend on the derivatives of the basic state. Furthermore, we assume that entropy has a 
maximum in thermodynamic equilibrium, that is a. = at equilibrium. According to 
the Second Law, the entropy function is negative semidefinite and increasing along the 
processes determined by Eq. @. Therefore, the entropy balance is given as 

(8) p s + V • J = a s > 0. 

Here p s — ps is the entropy density, J = j — ^ is the extra entropy flux and V- is the 
divergence operator. 
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Dissipation can be calculated according to Eq. (0 

n 

(9) T<r s = -A-g a = -Y,Ai9 l >0, 

i=l 

where the thermodynamic affinity A conjugated to the internal variable a is determined 
by Eq. (0 and • denotes a duality mapping onto the space of the internal variables as it is 
given by indexes. 

3.2. Fluxes and forces. We may recognize a simple force-flux structure in Eq. (O and 
identify A as a thermodynamic force and g Q as a thermodynamic flux. A general solution 
of this inequality can be given in the form of so called quasilinear conductivity equations 
(28l|23 

(10) g a {a) = -L(A(a),a)A(a). 

Substituting the latter into Eq. ^} results in a relaxation dynamics, because the conduc- 
tivity matrix L is positive semidefinite, as a consequence of Eq. (|9]). We should make a 
distinction of this quasilinear case, where the conductivity matrix is a function of the forces 
and the state space (L(A(ct), a)), and of the strictly linear case, where the conductivity 
matrix is constant (L = const). The strictly linear theory is always an approximation that 
can be valid only in the vicinity of equilibrium. 

The negative sign of the affinity is traditional and can be understood from the standard 
equilibrium formulation. 

3.3. Dissipation potentials. General dissipation potentials were introduced by Onsager 
Ell [30] in the case of strictly linear conductivity equations. From the point of view of the 
general quasilinear conductivity equations fllOt . the existence of dissipation potentials is 
connected to the validity of the general Gyarmati-Li reciprocity relations 13T1I321 . In our 
case, the Gyarmati-Li reciprocity relations in the so-called force representation [ 33 1 require 
that the derivative of the thermodynamic fluxes g Q with respect to the thermodynamic 
forces (affinities A), i.e., |^(A, a) should be symmetric. This means that there exists a 
dissipation potential Da{A, ex) with the property 

-^-(A,a) = g a {A,a). 

Formally, assuming a symmetric relation of the thermodynamic fluxes and forces, one can 
introduce the so-called flux representation 11331 supposing that the thermodynamic forces 
are functions of the thermodynamic fluxes. In this case, the condition of the existence of 
the dissipation potential is connected to the validity of the related Gyarmati-Li reciprocity 
relations, that is to the symmetry of the derivative ^^-{g a , a). As a result, we can get a 
flux related dissipation potential D g with the property 



'9 
dD, 



(g a ,a) = A(g a ,a) 



dg a 

In the strictly linear theory with a constant conductivity matrix L, these two representa- 
tions can be transformed into each other and both conditions coincide. Moreover, in this 
case the general Gyarmati-Li reciprocity relations become equivalent to the Onsagerian 
reciprocity relations. In the quasilinear case, they are not equivalent. It is also important 
to remark that the role of fluxes and forces is not completely interchangeable, because the 
fluxes are constitutive functions, and the forces are given functions of the constitutive state. 
Therefore, the validity of the flux representation is restricted mostly to the strictly linear 
theory. 
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3.4. Weakly non-local internal variables. We now release the local state hypothesis and 
assume that the entropy is a function of both the internal variables and their gradients 
(a, Va). Moreover, one can postulate ([3|, |8|) that the extra entropy flux is 

<"> J = " 

\dva / 

In this case, the entropy balance ([8]l together with the evolution equation (0 results in 

p s (a,Va) + V • J = 

dp s . dp s ( dp s 

-Va - V • -7—— 9a 



da <9Va \<9Va 

(12) = f^-V-|^W = -A. 9Q >0. 

\ da ova J 

Now the thermodynamic flux is the same as in the local case, but the thermodynamic force 
has been changed to the new weakly non-local affinity A = — jj^- + V • J^j- Therefore 
the evolution equation with quasilinear conductivity is 

(13) a = I Y^£_ V . 9Ps 



da dVa / 

It is remarkable that the right hand side of the evolution equation become nonlocal as well. 
The forces and the fluxes are compared in the following way: 

Local force: fg-(a), Nonlocal force: §^(a, Va) - V • -§^(a, Va), 
Local flux: g a {a), Nonlocal flux: g a (a, Va, V 2 a). 

Dissipation potentials can be generated as previously, but the previous local affinity A 
should be changed to the new weakly non-local affinity A. 

4. Internal degrees of freedom 

Evolution equations of other microstructural variables - the internal degrees of freedom 
- are generated by mechanical principles. It is assumed that the dynamics of an internal 
degree of freedom a is determined by a variational principle of Hamiltonian type with the 
Lagrangian C(a, a, Va) and, therefore, it is governed by a field equation of the canonical 
form (e.g. H): 

{ ' 5a ~ da dt \da J ' \d(Va) J ~ 

Here f Q is the dissipative force and is a functional derivative. 

The Lagrangian C is usually divided into a kinetic and a potential part as follows: 

(15) C(a,a, Va) = K{a) - V{a, Va). 

Substituting this relation into Eq. ( fT4l i, one can get the following particular form of the 
evolution equation for the internal degree of freedom a 

d 2 K .. dV dV 

da 1 da dva 

where we have rearranged the terms to get an apparent Newtonian form. 

One can consider higher order space derivatives and generate theories with higher or- 
der non-locality. On the other hand, one can apply the local state hypothesis, where the 
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Lagrangian does not depend on the space derivatives and the above equation simplifies. 
The natural boundary conditions for the internal degree of freedom are consequences of 
the variational formulation: 

(17) - ^ 



dVa gv 

where n is the normal of the boundary of the considered region dV at a given point and 
T is the surface "force" acting on the field a. The presence and the nature of natural 
boundary conditions is sometimes connected to the observability and controllability of the 
internal degree of freedom J5] . 

Let us observe that the above mechanical description of the dynamics of the internal de- 
gree of freedom postulates a Hamiltonian variational principle and immediately generates 
a second order differential equation in time ( [ToT i. However, there is a natural way to get two 
first order differential equations instead of a second order one introducing the Hamilton- 
ian through a suitable Legendre transformation. This transformation reveals some internal 
symmetries of the whole variational structure. 

In fact, introducing a generalized momentum 

dC 
Pa = ^-r, 

oa 

and a Hamiltonian H, defined by the partial Legendre transformation, 

p a a = C(a, a, Va) + H(p a , a, Va), 
we arrive at the first Hamiltonian equation 

(.8) d - " 

On the other hand, the field equation ([Pfl l can be transformed resulting in the second Hamil- 
tonian equation 

(19) Pq = -_ + V ' 



da \d(Va) 

Let us consider a particular form of the Hamiltonian equations in the case of the special 
Lagrangian given in Eq. Q5j. For simplicity, we also assume that the kinetic term K 
is quadratic, i.e. K(a) — ma /2, where m — const. In this case, the generalized 
momentum is p a = ma, and one can get <i(p Q ) = p a /m. The Hamiltonian follows as 

2 

H(a, p a ) — + V(a, Va). Finally, the corresponding special Hamilton equations are 
(20) a = Pa 



in 



dV ( dV \ 

4. 1 . Dissipation potentials. Some thermodynamic background can be added through the 
observation that the right hand side of Eq. (fT4l . the so-called dissipative force f Q , enters the 
entropy production because it generates dissipated power. The external mechanical power 
is included into the balance of total energy, then in the balance of internal energy (f2j), and, 
as a consequence, in the entropy balance (fjjl, as well (an other reasoning is based on the 
principle of virtual power (8] [34]). Finally, the entropy production can be written as 

Tcr s = f Q • a > 0. 
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The whole thermodynamic formalism of the previous section can be repeated on the ba- 
sis of the latter form of the entropy production, but now the thermodynamic fluxes are 
dissipative forces. However, here the thermodynamic formalism is added to the previous 
considerations of mechanical origin. For example, it is convenient to introduce dissipation 
potentials (with the above mentioned conditions) D f in the flux representation and put the 
dynamic equation ([Pil l into the following form [3 1 

5C _ 3D j 
5a df a 

5. Dual internal variables 

As we have seen, evolution equations for internal variables of state and for internal 
degrees of freedom are completely different (cf. Eqs. (0 and ([Toll). The same is clearly 
demonstrated recently on the example of continuum thermomechanics 0. 

The question arises whether it is possible to construct a single way of the derivation of 
evolution equations for both internal variables of state and internal degrees of freedom, or 
not? 

A first guess may be that the internal variables of state are special cases of the inter- 
nal degrees of freedom, because it looks like very easy to get a first-order time-derivative 
equation from the second-order one. However, this is not the case. On the contrary, the 
structure of Hamiltonian differential equations ( TT~8T > and ( fl9T i is very special. The essential 
part of the problem is that a second-order time-derivative equation can be generated by 
a Hamiltonian variational principle similar to the traditional variational principle in me- 
chanics without any further ado, but a first-order time-derivative equation only in special 
constrained cases (gyroscopic degeneracy (35)). This is a strong mathematical restriction, 
and any attempt to circumvent the problem has a price of loosing some parts of the nice 
Hamiltonian structure l20l [36ll . 

Another important observation is that the dissipative part of the dynamics is generated 
by dissipation potentials also in the second, mechanical method. It is remarkable that the 
mechanical generation of evolution exploits thermodynamical methods. 

In order to answer the formulated question, we introduce a dual internal variable and 
compare the arising evolution equations with the corresponding equations obtained in the 
previous sections. 

Let us consider a thermodynamic system where the state space is spanned by two scalar 
internal variables a, (3, Then the evolution of these variables is determined by the following 
differential equations 

(22) a = g a , 

(23) $ = g p , 

where the functions g a and g@ of the right hand side of the differential equations are con- 
stitutive functions and should be restricted by the Second Law of thermodynamics. The 
entropy inequality, the main ingredient of the Second Law, is the same as previously (Eq. 
dHJ). Let the domain of the constitutive functions (our constitutive space) is spanned by 
the state space variables and by their first and second gradients. Therefore, our constitutive 
functions g a ,g/3, J and p s are given as functions of the variables a, Va, V 2 a, (3, V/3, V 2 /3. 
This is a weakly nonlocal constitutive space with second order weak non-locality in both 
variables. 



INTERNAL VARIABLES AND DYNAMIC DEGREES OF FREEDOM 



9 



5.1. Liu procedure. The Second Law restricts the form of the possible evolution equa- 
tions (Eqs. (T22l i-(l23l>) and several exploitation methods of the Second Law can be imple- 
mented 11371 . Here we apply the procedure of Liu [38 1. 

The gradients of the constitutive functions g a and gp are constrained by the entropy 
inequality in the framework of a second order constitutive state space for both of our vari- 
ables as follows [39|: 

(24) Vd = V.g Q (a, Va, V 2 a,/3, V/3, V 2 /3), 

(25) V/3 = Vgp{a, Va, V 2 a,(3, V/3, V 2 /?). 

Let us introduce the Lagrange-Farkas multipliers A Q , Ay a , A/3, Ay/3 for Eqs. d22L d24"l i. 
(l23l and ( |25] |. respectively. Developing the partial derivatives of the constitutive func- 
tions, we use the the numbering of the variables of the constitutive state space spanned 
by (a, Va, V 2 a, /3, V/3, V 2 /3) as a convenient and short notation, e.g., d 3 p s — g^-! a ■ 
Therefore, Eq. ([8]) can be written as follows: 

p s + V • J - A Q (d - g a ) - A VQ (Vd - Vg a ) - A (/3 - g fj ) - A V/3 (V/3 - V 5/3 ) = 
<9ip s d + <9 2 p s Vd + d 3 p s V 2 a + <9 4/ 9 S( 3 + d 5 p s Vf3 + <9 6 p s V 2 /3+ 
6>i J • Va + <9 2 J : V 2 a + <9 3 J- : V 3 a + <9 4 J ■ V/3 + <9 5 J : V 2 /3 + <9 6 J- : V 3 /3- 

A a (d - </ a )- 

A Va (Vd - d iga • Va - 9 2 <7 Q : V 2 a - d 3 g a - : V 3 a - d i9a ■ Vf3d 5 g a : V 2 /3- 

d 6 g a - : V 3 /3)- 

A V /3(V/3 - <9i 3/3 ■ Va - : V 2 a - d 3 g p - : V 3 a - dtgp • V/3 - <9 5 <?/3 : V 2 /3- 

(26) 9 65/3 - : V 3 /3) > 0. 

Double and triple dots denote inner products of the corresponding tensors (scalars are 
formed). One can group together the coefficients by the derivatives of different orders 

a {dip s - A a ) +0 (d 4Ps - A p ) + 
Vd (d 2 p s - A Va ) + V/3 (d 5Ps - A V/3 ) + 
V 2 ad 3 p s + V 2 /3d 6 p s + 
V 3 a (d 3 3 + A Va d 3 g a + A vp d 3 g p ) + V 3 /3 (d 6 J + A Va d 6 g a + A V pd 6 g ) + 
V 2 a (d 2 3 + A Va d 2 g a + A V pd 2 gp) + V 2 /3 (d 5 3 + A Va d 5 g a + A V/ 3<9 5 ff/3) + 
Va (dxJ + Ava^iPa + A^pdigp) + V/3 (<9 4 J + A Va d 4 g a + A^pdigp) + 

(27) A a g a + A [jgp > 0. 
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We find the Liu equations as the coefficients of the derivatives that are not in the constitutive 



space 










(28) 


d 


dip s 


- A„ = 0, 




(29) 


$ 




- A = 0, 




(30) 


Vd 


2 p s 


- A Va = o, 




(31) 


V/3 


d 5 p s 


- A v/ 3 = 0, 




(32) 


V 2 d 


dsPs 


= 0, 




(33) 


V 2 /3 


dap s 


= o, 




(34) 


V 3 a 


5 3 J- 


f- A Va d 3 g a + A V /303ff/3 


= o, 


(35) 


V 3 /3 


d 6 J- 


f A^ a d 6 g a + A^pd 6 gp 


= 0. 



The residual dissipation inequality follows from Eq. (|27| | by taking into account Eqs. 



(01 J + A Wa dig a + A^pdigp) • Va + (0 4 J + A Va dig a + A V pd 4 gp) ■ V/3+ 
(0 2 J + A Va d 2 g a + Ay/3025/3) : V 2 a + (0 5 J + A Va d 5 g a + A v0 d 5 g p ) : V 2 /3+ 
(36) dip s g a + dip s gf3 > 0. 

It is easy to see, that higher than first order derivatives of the constraints do not give addi- 
tional restrictions, due to the simple structure of the evolution equations (l22l> - d23t . 

The solution of the Liu equations (f28T> - d33T> leads to the entropy density, which is inde- 
pendent of the second gradients p s (a, Va, V 2 a, /3, V/3, V 2 /3) = p s {ct, Va, (3, V/3), and 
its partial derivatives are the Lagrange-Farkas multipliers: 



A a = dip s = A Va = d 2 p s 



dp s 



A/3 = d 4 p s = A V/ 3 = d 5 p s 



dp s 



0Va' p 0/3' vp Jrs 0V/3 

Substituting the Lagrange-Farkas multipliers into Eqs. d34l i and d35l l, one can solve them 
and get the entropy flux 

J(a, Va, V 2 a, /3, V/3, V 2 (3) = ~d 2 p s g a - d 5 p s gp 



-Jo 



(37) 



dp s 
0Va 



9a 



dp s 
0V/3 



9f3 + Jo- 



Here the additional entropy flux term j is a function only of the variables (a, Va, /3, V/3), 
as well as p s . Substituting this expression into the dissipation inequality, we get finally: 



(38) 



<t s = {dip s - V • d 2 p s )g a + (d 4 p s - V • d 5 p s )gp + V • j = 

dp s 

0/3 



d£s_ _ dp s 
da dVa 



9o 



Y7 d ' Ps \ 

v -0v^ 90 



V jo > 0. 



Let us remember, that p s is independent of the second gradients of the variables, therefore 
there are no higher derivatives in Eq. d38l . Let us assume now that jo = 0. This is an 
assumption that we regularly consider in classical irreversible thermodynamics excluding 
the appearance of current multipliers and terms similar to that appeared in the Guyer- 
Krumhansl equation [40, 27]. We may recognize an Onsagerian force-flux system, where 
the thermodynamic fluxes and forces are 



a-force: A 
/3-force: B 



da 

OPs 

' OP 



v- 

dVa 
v dVP 



a-flux: 
/3-flux: 



9a 
9f3 
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Here we introduced a shortened notation A and B for the thermodynamic forces preserving 
the sign conventions of the previous sections. 



5.2. Evolution equations. Solution of the inequality (1381 1 can be represented as 

fc) £)(-t 

This corresponds to the inequality (l38l l if the Onsagerian coefficients Ii,l2,h2>hi are func- 
tions of the thermodynamic forces A and B and the state variables (a, Va, /?, V/3). Signs 
of the coefficients should be further restricted by the requirement of nonnegative entropy 
production. 

In what follows, we split the conductivity matrix into a symmetric and an antisymmetric 
part, introducing k and I instead of £12 and Z21 by 1 12 — I — k and I12 = I + k: 

(40) a = g a = kB- l t A - IB 

(41) $ = g = -kA-lA-l 2 B. 



The nonnegativity of the entropy production d38b results in the positive semidefiniteness of 
the symmetric part of the conductivity matrix L. Therefore, the following inequalities are 
valid for the coefficients 

(42) h > 0, h> 0, hh - I 2 > 0. 

5.3. Remark on dissipation potentials. We may introduce dissipation potentials for the 
dissipative part of the equations, if the condition of their existence is satisfied. However, 
in the considered case there is no need of this assumption, because here it is clear what be- 
longs to the dissipative part and what belongs to the nondissipative part of the the evolution 
equations. The terms with the symmetric conductivity contribute to the entropy production 
and the terms from the skew symmetric part do not. On the other hand, there is no need of 
potential construction, as we are not looking for a variational formulation. Moreover, the 
symmetry relations are not sufficient for the existence of dissipation potentials in general, 
as we have emphasized previously. In the case of constant coefficients (strict linearity), the 
dissipation potentials always exist for the dissipative (symmetric) part. 

5.4. Remark on the reciprocity relations. The reciprocity relations are the main results 
of the great idea of Lars Onsager connecting fluctuation theory to macroscopic thermody- 
namics [30l |4T1 l42l l43l . As it was written by Onsager himself on the validity of his 
result: "The restriction was stated: on a kinetic model, the thermodynamic variables must 
be algebraic sums of (a large number of) molecular variables, and must be even functions 
of those molecular variables which are odd functions of time (like molecular velocities)" 
ll42l . The Casimir reciprocity relations are based on microscopic fluctuations, too ll43l . 
We do not have such a microscopic background for most of internal variables. E. g., in the 
case of damage we think that the internal variables are reflecting a structural disorder on 
a mesoscopic scale. The relation between thermodynamic variables and the microscopic 
structure is hopelessly complicated. On the other hand, the Onsagerian reciprocity is based 
on time reversal properties of corresponding physical quantities either at the macro or at 
the micro level. Looking for the form of evolution equations without a microscopic model, 
we do not have any information on the time reversal properties of our physical quantities 
neither at the micro- nor at the macroscopic level. Therefore, we can conclude that lacking 
the conditions of the Onsagerian or Casimirian reciprocity gives no reasons to assume their 
validity in the internal variable theory. 
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Let us observe the correspondence of evolution equations for internal variables with the 
reciprocity relations by means of a few simple examples. 

5.5. Example 1: internal variables of state. Let us consider materials with zero anti- 
symmetric part of the conductivity matrix L (I = 0, k = 0). It is clear that the Onsagerian 
reciprocity relations are satisfied, and we return to the classical situation with fully uncou- 
pled internal variables: 

(43) a = g a = -l\A, 

(44) $ = g = -l 2 B. 

In this case the evolution equations for dual internal variables a and j3 are the same as in 
the case of single internal variable of state (Eq. ([T3~l>). 

5.6. Example 2: internal degrees of freedom. We now assume that all conductivity co- 
efficients are constant, and their values are l\ = I = 0, k = 1. This means that I12 = —hi, 
i.e., the Casimirian reciprocity relations are satisfied. 

For simplicity, we decompose the entropy density into two parts, which are dependent 
on different internal variables 

(45) p s (a, V«, P, V/3) = -K{J3) - W(a, Va). 

The negative signs are introduced taking into account the concavity of the entropy. Then 
the thermodynamic forces are represented as 

and Eqs. (|40b-(f4~Tb are simplified to 

(46) a = g a = B = K'{(3) 

dW dW 

(47) = 9l3 = -A + l 2 B = -— + V- — -hK'(j3). 

oa ova 

One may recognize that the obtained system of equation corresponds exactly to what was 
introduced in the case of dynamic degrees of freedom after transforming the equation of 
motion into a Hamiltonian form ( TT~8b and ( fT9l ). 

The transformation into a Lagrangian form is trivial if K has a quadratic form K(f3) = 
where m is a constant. Then we get exactly Eqs. (f20b - (f2TT> and the whole system 
corresponds to Eq. (5.14) in with the Lagrangian 

C(a, a, Va) = m— W(a, Va), and D(a, a) — — — at . 

Moreover, the entropy flux density (|37| | in the case of our special conditions can be written 

as 

(48) J = -|^) + jo. 

and one can infer that natural boundary conditions of the variational principle ( fTTT i corre- 
spond to the condition of vanishing entropy flow at the boundary, i.e., jo = 0. 

Therefore, the variational structure of internal degrees of freedom is recovered in the 
pure thermodynamic framework. The thermodynamic structure resulted in several sign 
restrictions of the coefficients, and the form of the entropy flux is also recovered. The 
natural boundary conditions correspond to a vanishing extra entropy flux. 
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5.7. Example 3: diffusive internal variables. Now we give an additional example to see 
clearly the reduction of evolution equations of internal degrees of freedom to evolution 
equations for internal variables and the extension of the later to the previous one. 

We keep the values of conductivity coefficients (i.e., 1% = I = 0, k = 1), but assume 
that both K and W are quadratic functions 

where a and b are positive constants according to the concavity requirement. 
In this case, the evolution equations d46l i. ( f4Tb reduce to 

(49) a = g a = B = K'((3)=bp 1 

dW dW 

(50) (3 = 9l3 = -— + V- — ~l2K , (f3)=aAa-l 2 b[3. 

oa ova 

Substituting /3 from Eq. d49l into Eq. (f50]l, we have 

1 li 

(51) — a+— d = Aa. 

ab a 

which is a Cattaneo-Vernotte type hyperbolic equation (telegraph equation) for the internal 
variable a. This can be considered as an extension of a diffusion equation by an inertial 
term or as and extension of a damped Newtonian equations (without forces) by a diffusion 
term. 



6. Discussion 

In the framework of the thermodynamic theory with dual weakly nonlocal internal vari- 
ables presented in the paper, we are able to recover the evolution equations for internal 
degrees of freedom. 

We have seen that the form of evolution equations depends on the mutual interrelations 
between the two internal variables. In the special case of internal degrees of freedom, the 
evolution of one variable is driven by the second one, and vice versa. This can be viewed as 
a duality between the two internal variables. In the case of pure internal variables of state, 
this duality is replaced by self-driven evolution for each internal variable. The general case 
includes all intermediate situations. 

It is generally accepted that internal variables are "measurable but not controllable" (see 
e.g. 11251 ). Controllability can be achieved by boundary conditions or fields directly acting 
on the physical quantities. We have seen how natural boundary conditions arise considering 
nonlocality of the interactions through weakly nonlocal constitutive state spaces. 

As we wanted to focus on generic inertial effects, our treatment is simplified from sev- 
eral points of view. Vectorial and tensorial internal variables were not considered and 
the couplings to traditional continuum fields result in degeneracies and more complicated 
situations than in our simple examples. 

A number of systems with these degeneracies are called gyroscopic. One can distinguish 
a least two types of them. In thermodynamics this is a customary nomination for skew 
symmetric Casimir type couplings of thermodynamic interactions, because in this case the 
entropy production is zero. For example a single second order tensor internal variable £ 
can couple to mechanical interactions for isotropic materials, too. In this case the relevant 
part of the entropy production from Eq. (O can be written as: 



(52) 



Ta s = (t - t rev ) : Vv - H : | > 0. 
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Here S = — is the related affinity. If £ is a variable that changes sign under time 
reversal then there is a pure skew symmetric coupling between the two variables and the 
entropy production is zero. This is a valid and important example in rheology (see |6| 
for second order traceless tensors as internal variables). The above example is sometimes 
called as gyroscopic (e.g. [44 J), because the corresponding thermodynamic forces and 
fluxes are orthogonal. 

On the other hand the gyroscopic systems in mechanics are typically characterized by a 
Lagrangian with the following scalar product term: 

L(q,q) = a(q) • q, 

where a is a vector field, that depends on the generalized coordinates q of the mechani- 
cal system and has the same dimension as coordinate space. The dot denotes the scalar 
product. The corresponding Euler-Lagrange equation follows: 

L>Aaq = 0. 

Here the derivative D A a = Da — (Da) T . This system is degenerate, because the corre- 
sponding Legendre transformation is not invertible, the related Hamiltonian is identically 
zero. 

These two examples show well that skew symmetric couplings are not always related 
directly to inertial effects and indicate two directions, one into mechanics and one into 
thermodynamics, where our method can be generalized. Let us mention here the related 
pioneering works of Verhas, where skew symmetric conductivity equations appear in dif- 
ferent inspiring contexts [45, 46). 

Finally, let us mention that the idea of constructing a unified theoretical frame for re- 
versible and irreversible dynamics has a long tradition. The corresponding research was 
not restricted to the case of internal variables and was looking for a classical Hamiltonian 
or a generalized variational principle that would be valid for both dissipative and nondissi- 
pative evolution equations (see, e.g., (33]|47]|48l[36l and the references therein). 
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